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I. INTRODUCTION 
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We study muonium-antimuonium conversion and fj, + e~ — * n~e + scattering within two different 
lepton-flavor-violating models with heavy neutrinos: model I is a typical seesaw that violates lepton 
number as well as flavor; model II has a neutrino mass texture where lepton number is conserved. 
We look for the largest possible amplitudes of these processes that are consistent with current 
bounds. We find that model I has very limited chance of providing an observable signal, except if 
a finely tuned condition in parameter space occurs. Model II, on the other hand, requires no fine 
, , tuning and could cause larger effects. However, the maximum amplitude provided by this model is 

still two orders of magnitude below the sensitivity of current experiments: one predicts an effective 
coupling G MA j up to 10~ 4 Gf for heavy neutrino masses near 10 TeV. We have also clarified some 
, discrepancies in previous literature on this subject. 

> 

Muonium is a hydrogcnlikc system of an electron (e ) bound to a positive muon (/i ) and antimuonium is the 
bound state of the respective antiparticles (fi~e + ). It is formed as muons stop inside matter, although the fraction 
of stopped muons that actually form muonium depends very much on the material itself. Since it has no hadronic 
constituents, muonium is an ideal place to test electroweak interactions. In particular, since conversion of muonium 
into antimuonium is completely forbidden within the standard model (SM) because it violates lepton flavor, its 
,H« observation will be a clear signal of physics beyond the SM. This phenomenon was first suggested almost 50 years 
Oh, ago [| and since then several studies 0, Q and experimental searches £| have been done. Many extensions of the 
SM could cause muonium to antimuonium conversion [5j (left-right symmetric gauge models, extra neutrinos, extra 
D | higgses, SUSY, bileptonic gauge bosons, etc). 

• Here we want to study the possibility that heavy Majorana neutrinos in models of seesaw type could cause an 
observable signal of muonium-antimuonium conversion. In previous works we have studied the prospects for these 
models to generate tau decays that violate lepton flavor [|| . The main feature of these models is that they only extend 
the standard model by the existence of extra neutral leptons. The most natural scenario is then that lepton families 
will mix (thus destroying lepton flavor conservation) and that neutrinos will acquire masses. In order to insure that 
the standard neutrinos remain light, the texture of the neutrino mass matrix should favor a seesaw type of mechanism. 
In these models, muonium to antimuonium conversion appears at one loop level in the perturbative expansion (box 
diagrams -see Figs. 1 and 2) and it is thus a very small effect. 

A simple estimate of this conversion probability within the models in question shows that it is too small to be 
observed in foreseeable experiments 0j. Another work along these lines making a rough estimate shows that 
present experimental bounds set a lower bound on large Majorana masses at around 10 1 TeV. However, we find 
that the validity of the calculation (perturbation theory to one loop) breaks down for masses above that same scale. 
Therefore one can make clear predictions for masses up to that scale only. We analyze the process in more detail within 
the one-loop approximation, performing a comprehensive search for the bounds in the currently allowed parameter 
space. We confirm the conclusion that seesaw models give a conversion probability that is still too small to be at the 
reach of experiments to date. We also find that typical seesaw models of heavy neutrinos are in most cases totally out 
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of reach of any foreseeable experiment, except in very special, fine-tuned regions of the parameter space. In contrast, 
models with heavy neutrinos that conserve lepton number give conversion probabilities closer to present experimental 
sensitivities for natural values of parameters. 

In section II we give the generalities of muonium-antimuonium conversion; in section III we calculate the process 
within the models in question and do a comprehensive search of the presently allowed parameter space, looking for the 
maximal possible conversion probabilities. In section IV we give brief account of the unbound collision ji + e~ — ► /i~e + 
and its prospects. In section V we state our conclusions. Details of the seesaw models and loop calculations are 
included as appendices. 



II. MUONIUM TO ANTIMUONIUM CONVERSION 



Muonium \M) is a nonrelativistic Coulombic bound state of fi + and e~ , and antimuonium \M) a similar bound 
state of n~ and e + . The nontrivial mixing between \M) and \M) comes in our context from the non- vanishing LFV 
transition amplitude for e~ fi + — ► e + /i~. In a model with heavy neutrinos, this amplitude appears at one loop (box 
diagrams) in the electroweak interaction extended with extra (heavy) neutrinos. This transition is usually expressed 
in terms of a local effective Hamiltonian density with an effective coupling G MM '- 

H cS {x) = ^ [7l(x) 7 "(l- 75 ) e (x)] [p(x) 7a (l-y 5 )e(x)] . (1) 

Due to this interaction, \M) and M) are not mass eigenstates. In this basis, the mass matrix has non-diagonal 
components: 



'■MM 



(M\ / d 3 x H(x)|Af), (2) 



where the (non-relativistic) bound states are chosen with unit norm in a volume V, namely (M(P)| Af(P')) = 
(2n) 3 S 3 (P — P')/V. If we express the bound state in terms of its constituents with a momentum distribution f(p): 

\M(0)) = J ^f(pK e)tb -p\°)> (3) 
and neglect the momentum dependence in the spinors, the mass mixing element is: 



l MM — 10 



V2 



d 3 p 
(2^)3 



u^)Ll a u {e)L v [ll)L -f a v {e)L (4) 



Since angular momentum is conserved, the only non vanishing elements are those where the spin of the initial and final 
bound state are the same. For both the singlet or each of the triplet components, the spinor product has the same 
value 2m^m e 9]. Additionally, the integral of f{p) is the spatial wavefunction at zero distance, | J d 3 p/{2it) 3 f(p)\ 2 = 
| , 5(0)| 2 /(2m Al 2m e ). The mass mixing element thus reduces to 

m _ w g M m 1*(Q)| 2 _ l( ,G MM (m? f5l 

m MM -W— — (5 ) 

where /i rs m e is the muonium reduced mass and a the fine structure constant. In vacuum, the diagonal mass terms 
are equal, tumm = m MM- The mass eigenstates are then the simple combinations \M\) — (|M) + |M))/V2 and 
\M 2 ) = (\M) — \M))/y/2 and the mass eigenvalues differ in Am = mi — m 2 = 2 m M ^. This is the typical case of a 
quantum system of two states that oscillate into each other under time evolution, like the K° — K° system in hadrons. 
However, unlike the latter, in practical terms muonia do not have time to show an oscillation, because their oscillation 
time is much longer than their decay time. So one can only hope to observe the mixing phenomenon by measuring 
the probability that a state that starts as a Muonium (fi + e~) decays as an Antimuonium (decay into a high energy 
electron and low energy positron). So, the state that starts as \M) at t — (call it \M(t))) evolves as: 

\M(t)) = |Mi)(M 1 |M)e- i mit + \M 2 ){M 2 \M)e- 1 ™ 2t (6) 

and will decay as a M with a total probability: 
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P(M -> M) 


Gmm /Gf 


Experiment 


< 2.1 x 10" 6 


< 0.29 


Huber et al. (1990) 


< 6.5 x 1(T 7 


< 0.16 


Matthias et al. (1991) 


< 8.0 x 10~ 9 


< 1.8 x 10~ 2 


Abela et al. (1996) 


< 8.3 x 10" 11 


< 3 x 10" 3 


Willmann et al. (1999) 



TABLE I: Experimental results (see Ref. Q) for the conversion probability and the deduced upper bound for the coupling 
Gmm (this bound on the last row is weaker than what Eq. JSJ gives, because of magnetic field suppression). 



where r is the muon lifetime. Eq. Q does not take into account the effect of static electromagnetic fields in materials, 
which break the degeneracy tumm = m MM an d cause an extra suppression on the probability 0, an important effect 
in some experiments. In any case, we see from Eq. Q that the conversion probability is in general very small and 
proportional to \G M m\ 2 : 

P(M -> M) = 64 G F^ m " t2 (^0-) 2 = 2.64 x 1CT 5 {^f) 2 ■ (8) 

The experiments so far have gradually been improving on an upper bound for P(M — > M) and thus for G M! ^ (See 
Table 1), but no positive signal of this lepton flavor violating process has ever been observed, which is consistent with 
the null prediction of the standard model. All new physics in this process enters in the parameter GjvfM- We will 
now examine the prospects for G M ^ in extensions of the standard model with extra heavy neutrinos originating in 
seesaw mechanisms. 



III. PREDICTION OF THE SEESAW MODELS 



We want to explore the prediction of lepton flavor violating seesaw models for muonium-antimuonium conversion. 
The details of the models are given in Appendix A. The main features of these models are that (i) they contain 
heavy neutrinos, (ii) in their mass basis, the heavy neutrinos couple to charged leptons of different generations under 
weak interactions (just like the Cabbibo-Kobayashi-Maskawa prescription in the quark sector), thus inducing flavor 
changing processes (iii) in general the neutrino masses are of Majorana type. 

The mixing matrix Bn, which is defined in Appendix A, is the coefficient of the left-handed current that connects 
a charged lepton of flavor I = e, fJ-, ■ ■ ■ with a neutrino of flavor i = 1, 2, . . . In a model with Nl charged leptons, 
the first Nl neutrinos are light (standard) and the rest are heavy. Here we will only consider Nl = 2 to make the 
estimates. Of course, general lepton flavor mixing would induce many kinds of processes that are not observed in 
the real world. The non-observation of those processes impose severe bounds on the values of Bu. Here we use the 
results of the analysis of Ref. (01), which establish bounds on combinations of mixing elements characterized by 
(s^ ! ) 2 = J2i=h \ B n\ 2 > where h indicates heavy neutrinos, singlets under SU(2)l- 

(s^) 2 < 0.005, (s^) 2 < 0.002, (s^) 2 < 0.010. (9) 

As we said, in these models, muonium-antimuonium conversion is induced by the process fi + e~ — > fi~e + at one-loop 
via box diagrams (see Figs. 1 and 2), where the neutrinos appear as internal quanta. The effective coupling G^m 
is a function of the internal neutrino masses m, and mixing parameters Bu that arises from the integrals of the box 
diagrams: 

g mm/ g f = -^- 51 [ 2B ti B tj B ^ B M FBox{xi,Xj) - (B* ei ) 2 (B^ j ) 2 X*X j G Box (xi,Xj) , (10) 

i,j—all v' s 

where aw — ff 2 /47r is the weak SU(2)l fine structure constant, Xi = m 2 /M w are the squares of neutrino masses inside 
the loop in units of the W boson mass, Xi are "creation phases" that appear in Majorana fields 0], and FboxC^ij Xj) 
and G-Box{xi,Xj) are the loop functions that appear integrating the "Dirac boxes" of Fig. 1 and the "Majorana boxes" 
of Fig. 2, respectively (external masses and momenta inside the integrals are neglected as compared to loop momenta). 
These functions, symmetric under the interchange Xi < y Xj , are: 

if/ XiXj\ ( 1 \ / 1 Xi 

FBox{Xi,Zj) = < 1H — h — r~mXi \-2xiXj h — r^-ma;, 

Xi-Xj{\ 4 /VI -Xi (l-Xi) 2 J ''V 1 -^ {l-Xi) 2 
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(4 + XiXj) 



-(i-j) 
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(12) 
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FIG. 1: "Dirac boxes": box diagrams for generic (Dirac or Majorana) neutrinos (each wavy line is either a W or a charged 
Goldstone boson if we work in the Feynman gauge). 
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FIG. 2: "Majorana boxes": box diagrams for Majorana neutrinos only (each wavy line is either a W or a charged Goldstone 
boson). 



We should point out that we have a discrepancy in the global sign in front of Gbox in Eq. IjlOfl compared to Ref. 
|l2j |. It is important to remark that this global sign can not be changed by simple redefinition of fermion phases. For 
example, one can redefine the phases of the Majorana fields and choose \ = +1 for all i, or equivalently one can 
absorb those phases in the mixing matrix Bu, but it is impossible to choose the phases so that A*Aj = —1 for all 
It is also important to notice a discrepancy with Ref. (0), where the authors claim that the contribution from the 
two diagrams of Fig. 2 with Majorana neutrinos (i.e. the Gbox term) cancel because of Fermi statistics. Actually, 
when one carries through the calculation, the minus sign from Fermi statistics is compensated by the inner structure 
of the diagrams and at the end both diagrams are equal and they add! Since Gbox dominates over -Fbox for large 
neutrino masses, this could be an important issue in the estimates. To clarify these matters, we include the derivation 
of Eq. (fTO)) in Appendix B. 

Using the unitarity of the mixing matrix, Eq. (|10|) can be rewritten as a double sum over heavy neutrino flavors 
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only 

N L +N H 

Gmm/Gf = ^f- V {2B* eI B* eJ B^B^j [F BaK ( Xl , xj) - 2F Box (^, 0) + F Box (0, 0)] 

OZ7T ' — ' l 

I,J=N L +1 

-{B^fiB^jfXjXj G Bok ( Xi ,xj)} , (13) 

where light neutrino masses are neglected. 

Table I shows that in order for muonium to antimuonium conversion to be observable in at least the most sensitive 
experiment to date, the coupling G mm /Gf should be at least 0(1CP 3 ). A rough estimate for the largest value 
expected for G MM from seesaw models of heavy neutrinos, and consistent with the present bounds of Eq. JJ]J, can be 
obtained considering that the functions -Fbox and Gb ox of Eqs. (fTTJ) and (fT2")l are growing functions for large neutrino 
masses and that the mixings in seesaw models are in general decreasing functions of the neutrino masses. For neutrino 
masses of a few TeV, -Fbox and Gbox are of order 10 1 and 10 3 respectively and we will assume that for those masses 
the mixings can still be near the maxima shown in Eq. ©. In that case the process is dominated by Gbox and we 
obtain 

Gmm/Gf ~ ^ G Bo x x ( S ^) 2 (.s£) 2 ~ 1(T 6 (14) 

which is 3 orders of magnitude smaller than present observable values. 

If we want to search for larger possible values for G MMl we should assume larger masses for the heavy neutrinos. 
However, the assumption of larger masses may enter into conflict with the choice of mixing parameters at their current 
upper bounds [Eq. |J3J]. As mentioned earlier, in seesaw models the mixings decrease as the neutrino masses are taken 
to be larger. If the neutrino masses are very large, the mixings will be forced to be considerably smaller than their 
current upper bound, and so assuming larger masses may not help to get potentially larger conversion probabilities. 
This effect was not an issue in Ref. || because they assumed rather large Dirac lepton masses (m^ ~ My/)- With this 
assumption one naturally predicts rather large conversion probabilities, however enters into conflict somewhere else: 
the prediction for the light neutrino masses, m u ~ rri^/Mh ~ 10 2 MeV, now turns out to be too large. Therefore, 
one needs to refine these estimates to be consistent with all present bounds. Here we do it by finding the relation 
between the neutrino masses and the mixings, following the whole process of mass matrix diagonalization from the 
start within each model. Keeping these consistency conditions, we search for the largest possible values for G M j^ in 
the scenarios of models I and II. 



A. Maximal conversion amplitude in Model I 



The key problem is to find the neutrino mass matrix M. in Eq. (|A5|) and the mixing parameters Bn in Eq. i|A8(l . 
constrained by Eq. @, that maximize the effective coupling G MM . 

The general form of M. for two families in Model I can be described as: 



M = 



rriFi 



mn = 




m M 




(15) 



where a, 6, c, d are real Dirac masses, 5j (J = 1,2) are two independent physical phases, and Mi and M2 (with 
Mi < M2) are large masses. The mixing parameters Bu are related to the transformation matrix U of Eq. I|A6() . 
which can be written approximately as a product of a seesaw block matrix U s and a light-sector mixing matrix V' 



u = V^U* 



U s = 



(hx2 - \m D mJrrV D ) 



V? 
7 2 x: 



cost* smtf e 
sin e le cos 8 



(16) 
(17) 



where e is a CP phase. The transformation 1A6(I with U s alone gives an approximately block-diagonal mass matrix if 
we neglect the correction terms 0{m 2 D /m 2 M ) and 0{rr? D jm} M \. 



U S MU 



-m D m^m T D (l + 0(m 2 D /m 2 M )) 0(m 3 D /m 2 M ) 



0(m 3 D /m 2 M ) 



m M (1 + 0(m 2 D /m 2 M )) 



(18) 
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The matrix , Eq. (|17|l . then generates a congruent diagonalization of the light-sector matrix m^m^m^. In , 
8 = and 7r/4 correspond to zero and maximal v e -v^ L mixing. According to solar neutrino experiments, we have 
maximal v e -v {i mixing, 6 = 7r/4. If we demand that this maximal mixing be obtained independently of the large 
masses Mi and M 2 , then straightforward algebra gives the following simple restrictions on the light Dirac sector 
parameters: 



2 2 
a = c , 



b 2 =d 2 , S 1 = S 2 (=5). 
The value of 5 can be restricted to —tt/2 < 8 < tt/2. The two light neutrino eigenmasses are then 



Mi 



Ah 



1I/2 



+ +2 cos(2<5) 



Mi M 2 



± 



ac bd 
M~i + M~ 2 



and the heavy-to-light mixing parameters (s^) 2 = Ylh=3 \Bih\ 2 and B* h = Um are 



(19) 



(20) 



(^) 2 = (*ln 2 = j~2 + J~2 (= s i < 0-002) 



5. 



//1 



-1 J "2 

a/Mi , (6/M 2 )' 
(c/Mi)e-« , d/M 



- i j 



£ = e, n ; ft, = 3, 4 



(21) 
(22) 



The elements of the diagonal matrix A of Eq. (| A6|) generally are complex phase factors; here we have A3 = A4 = +1 
while Ai, A 2 ^ +1. 

In order to reduce the parameter space, for every given Mi and M 2 we choose the values of the Dirac sector 
parameters (a, 6, c, d, Si, 8 2 ) that maximize the effective coupling |G Mi ^| [Eq- Q13[l]. and then explore the Mi,M 2 
parameter space on top of that condition. On the basis of Eq. I|13(l . we can set an upper bound on |G mj q-|, for given 
Mi and M 2 : 



oZtt 



\I=3 



7=3 



J=3 



(23) 



- 32tt V B 



G F 



7=3 



,7=3 



OZ7T 



(24) 



where -F Bo x = 2 max|F B ox(£7,£j)-2F B ox(2;7,0)-l-i 7 B ox(0, 0)|, and G Bo x = max|G Bo x(a;7, Xj)\. The maximum of G MJ tj 
is reached as it approaches this upper bound and simultaneously the maximal values (s^f ) 2 = (s^) 2 = ( s z,)max = 
0.002, Eqs. © and 0J, are reached. 

Now, there are two distinct cases into which the maximum mixing condition i|19|) can be divided: 

• Case (a): c = ±a and d = ±b (and 5 X = 8 2 = 6). In this case, m„ M > (a 2 /Mi + 6 2 /M 2 ) by Eq. 
and by Eq. (|21|) consequently s\ < m v /Mi < 3 ■ 10 -11 , where the latter inequality was obtained by ta king 



Mi > 100 GeV and m„ < 3 eV as suggested by present neutrino experiments (m 
\m 2 — m 



< 3 eV, Ref. 2 



'soi 



lO" 10 cV 2 0, Refs. |23|). This value of s| extremely suppresses G M jg, as seen by 
Eq. 124|) . so we will disregard this case for our purposes. 

Case (b): c = ±a and d = =F& (and <5i = <$2 = S). In this case, by Eq. H20(l 



?7?.,, 



a 

Wi 



M 2 



-2 cos(2<5) 

Mi M 2 V 7 



1/2 



± 



a 

Mi 



M 2 



(25) 



',„, - -m 2 c > 4|a 2 /Mi -fe 2 /M 2 | 2 , i.e., |a 2 /Mi - 6 2 /M 2 | < (1/2) (Am^) 1 ^ „ i -8 e V, 
which is very small compared to m e ,m ll . This means that solar neutrino experiments imply an almost exact 
relation: 



a 2 /Mi = b 2 /M 2 , 



(26) 
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which in turn implies the following expression for the mass: 

m v jM 1 = (a 2 /M 1 2 )V2(l + cos25). (27) 

This is a very small quantity, of order 10" 11 . Here we see the typical behavior of seesaw models: since m v /Mi ~ 
a 2 /M 2 and the mixing of heavy-to-light sectors is Bih ~ a/M\, the latter vanishes as the disparity between 
light vs. heavy neutrino masses intensifies, namely m„ <C M\. 

However, we do not want a/Mi to be too small, because it would suppress the mixing parameter s\ [see Eq. 
1)21(1] . Therefore, in view of Eq. (|27|l we set 6 — > 7r/2, which keeps the mass disparity without the vanishing 
of mixing. This is a rather finely tuned condition, but we admit it because we are only looking for maximal 
transition probabilities. Then, using (|21|l . (|26|) and the defining condition of this case, namely c — a and d = —b, 
we have set all the parameters of the Dirac sector for maximal conversion probability in terms of just M\ and 
M 2 : 

a = c = A/i(s L ) max /Vl + Mi/M 2 , b=-d = a^M 2 /M l , S — > it/2. (28) 

We then proceed to scan the values of |G MJ ij| as a function of Mi and M 2 restricted to the conditions of Eq. I|28|l . 
We find, interesting enough, that the Majorana boxes, i.e. terms proportional to Gb x(xi,xj), do not dominate the 
amplitude for large masses -at least not under condition (|28|l - even though the functions Gb x(^j, %j) are individually 
larger than Fb Q x(xi, xj). Instead, the mixing elements that multiply these functions are such that both Fb ox and 
Gbox terms are of the same order and interfere. Indeed, if one expands the functions Fbox and Gbox in the amplitude 
(f 1 31) for asymptotic values of xi and xj, the leading terms come naively from Gbox; however, these are multiplied 
by a factor proportional to (to„ /Mi) 2 arising from the mixing elements, which totally suppress them, leaving the 
leading terms of Fbox and the subleading terms of Gbox at the same level of interference. 

Fig. 3 shows the graphs for the maximal coupling | Caj'jv/ I as a function of Mi and M 2 , which is consistent with all 
present constraints within model I. 




2~ 4 6 8 ^10 

M2/M1 



FIG. 3: Maximal effective coupling G M m/Gf for muonium to antimuonium conversion in Model I, as a function of the ratio of 
Majorana masses M2/M1, for different masses Mi: from the lower to the upper curve, Mi = 100 GeV, 300 GeV, 1 TeV, 3 TeV, 
10 TeV. The dips reflect an accidental cancellation of the leading terms of Dirac and Majorana box diagrams for particular 
values of Mi and M2. 



We see that the amplitude is indeed a growing function of Mi , M 2 (we also see that for given Mi , local maxima 
are reached if M 2 « 2Mi). We do not push Mi,M 2 higher than around 10 TeV, because the perturbative expansion 
to one loop breaks down above that scale, as the authors of Ref. 01 point out. They stated that the perturbative 
expansion is valid provided M 2 J^e \Beh\ 2 < ^M^/aw- In our case, this condition becomes: 

M 2 

MiM 2 < ~ (10 TeV) 2 . (29) 

a w s L 
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From Fig. 3, we see that the maximum expected value of G M m/Gf is around 10~ 5 — > IO4, which is still almost two 
orders of magnitude below the level of sensitivity of current experiments. To close the analysis, we may comment on 
the finely tuned parameter <5 — > tt/2. It is clear that if we do not impose this condition, then the mixing elements s\ 
are ~ a 2 /Mf ~ m v /Mi. For m Vtl ~ 190 KeV, which is the upper bound and M x ~ 10 TeV, we get s\ ~ 10~ 8 , 
which signifies a suppression of 6 orders of magnitude in the amplitude with respect to the maxima found in our 
analysis. 



B. Maximal conversion amplitude in Model II 

In Model II, Eq. I|A12|) . with CP-conserving mass terms, the submatrices are real and can be taken in the form 



m D 



M 







m D 



mju 




m M 




Again here Mi and M2 (Mi < M2) are of the order of the heavy Majorana masses, and a,b,c,d ~ Mj\ s L 
diagonalization of matrix M. [see Eq. l|A6|l , but this time with U orthogonal] we get the neutrino masses 

mi = m 2 = , 

1 



(30) 
■ By 



7713 = 7714 = Mi 
m 5 = me = M 2 



1 + r 



1 a 2 + c 

2 M? 
lb 2 ' <* 
2 



M| 



(31) 



where the dots represent terms ~ 1/Mj. The phase factors in Eq. 1|A6[) are A3 = A5 = — 1, Ai = A2 = A4 = \q = +1 
and the mixing elements B^ arc real, fulfill the symmetry i?; 3 = B14 and Bi 5 — Bi 6 (for I — e,fi), and are given 
explicitly by: 



Bih = 



V2 



-a/Mi , -a/Mi , b/M 2 , b/M 2 
-c/M x , -c/Mi , d/M 2 , d/M 2 



O (m 3 D /m 3 M ) 



= e, )jl ; h = 3, 4, 5, 6 



(32) 



Notice that the two "light" neutrinos are actually massless, although they can easily acquire small masses by 
introducing small nonzero entries in the lower right 2x2 block of the mass matrix M.. The two degenerate pairs of 
heavy Majorana neutrinos can be interpreted as two heavy Dirac neutrinos. The symmetry in the masses, mixings 
and phase factors imply that the contribution to G MA j from the Majorana boxes of Fig. 2, [i.e. the terms with 
Gbox(xi, xj) in Eq.^5], cancel out. This contrasts with Model I where such terms are important, and it is related to 
the fact that Model II conserves total lepton number while Model I does not. 

As before, we search for the maximum of \G M m/Gf\- Since the Gbox(%i, xj) terms cancel, and the mixing elements 
Bih are real in absence of CP violation, we can write 

6 

\G M m/Gf\ = S^xl V \2B eI B eJ B llI B ll j [F Box { Xl , xj) - 2F Box { Xl , 0) + F Box (0, 0)] 

OZ7T I ' — ' I 

I,J=3 

Fbox ( / , B e iB^i 



< 



aw 
32tt 

aw 



< ^Box(4 C ) 2 (^) 2 



(33) 



where Fbox = 2max|i 7 Box(^/, xj) — 2Fbox(^/,0) + Fbox(0, 0)|. The last inequality is Schwarz's inequality, which 
saturates when the following proportionalities are reached: 



( ?3 



jU3 



e4 



^4 



(5 



fi5 



B 
B 



< 6 



//G 



(34) 



The first and the third equalities are fulfilled automatically. The second equality, however, can be shown to be fulfilled 
only when ad = be. Thus, the maximum of IGj^m] is reached when (s v £) 2 and (s^ 1 ) 2 are at their upper bounds, 
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Eqs. © 



Ml 



Ml 



Ml 



Ml 



0(m%/m 4 M ) 



0.005, 
0.002, 



(35) 
(36) 



and the condition ad = be is simultaneously fulfilled. This implies, for given values of high mass parameters Mi and 
Ma, three relations for the four low mass parameters a, b c and d. Thus, we still have one degree of freedom for 
fixing the values of low mass parameters, without much effect on the conversion probability. The simplest additional 
relation seems to be the symmetry assumption b = c. All this then gives the optimal choice of low mass parameters 
a, b, c and d in Model II: 



M, 



b = c = a X (#/#), <i = ax( S ^f, 



(37) 
(38) 



where s^ e and Sj? are taken at their respective upper bound. Again, we have reduced the parameter space to just 



two quantities, Mi and M^. Fig. 4 shows the graphs for the maximal effective coupling | a/aj I 
and M2, which is consistent with all present constraints within model II. 



as a function of M\ 



le-05 



le-06 



le-07 



le-08 



6 

M2/M1 
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FIG. 4: Maximal effective coupling G M m/Gf for muonium to antimuonium conversion in Model II, as a function of the ratio 
of Majorana masses M2/M1, for different masses Mi: from the lower to the upper curve, Mi — 100 GeV, 300 GeV, 1 TeV, 3 
TeV, 10 TeV. 



We see that in this model the amplitude G M ^ is an increasing function of Mi, but is rather insensitive to Mi 
{Mi > Mi by definition). Comparing Figs. 3 and 4 we see that in model II the maximum probability can reach 
slightly larger values than in model I. However, there is a stronger distinction between the models: here the maxima 
are reached for natural values of the parameters, while the maxima in model I are only reached for a finely tuned 
condition, which circumvents the seesaw suppression of the mixing; away from this finely tuned condition, the mixing 
suppression is very strong, pushing the amplitude of model I several orders of magnitude below the maximum. 
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IV. PROSPECTS FOR FREE e ~ fj.+ -» n~e + COLLISIONS 

As an aside from this calculation, one can notice that the same amplitude obtained from the box diagrams can be 
used to determine the cross section for free e~/i + — > /i~e + collisions. The result is simply: 



a(e n + -> n e+) = y-^HG MA7 /G 



4tt a 2 „ ,2 



|GmmI/Gf 



2 



~ 3x10 teA^^J fbanL (39) 

This expression is applicable for CMS energies well above but below Mw, i.e., for 1 GeV < -Js < 10 2 GeV. 
In this kinematical regime, the cross section is clearly too small to be observed in the foreseeable future, especially 
considering that both expressions in parentheses are at most of order unity. For this reason, we will comment on this 
process no further. 



V. CONCLUSIONS 



We have looked for the maximum probabilities of muonium to antimuonium conversion that can be reached in 
models where the required lepton flavor violation is caused by heavy neutrinos, keeping consistency with all present 
bounds. The estimates are done within two alternative models. In both models the amplitudes are in general growing 
functions of the heavy neutrino masses. Our estimates of the amplitudes are thus limited by two effects: (i) suppression 
of flavor mixing matrix elements as the masses get larger and (ii) the breakdown of perturbative expansion for masses 
above a few tens of TeV. 

Model I includes extra right handed neutrino singlets and constitutes a typical seesaw model for the neutrino masses 
and mixings, which violates lepton flavor as well as lepton number. Model II includes extra right handed and left 
handed neutrino singlets in such a way that lepton flavor is violated, yet lepton number is conserved. We found that in 
Model I the Majorana boxes give the dominant contribution over the Dirac boxes in the most typical scenarios, yet the 
maximum probabilities are obtained in finely tuned regions of the parameter space, where Dirac and Majorana boxes 
are comparable and even tend to interfere destructively. As a consequence of this effect, the maximal conversion rates 
are not as large as one could have naively estimated. The maximal couplings thus obtained are G M m/Gf ~ 10 -5 , 
which means conversion probabilities 4 orders of magnitude below present experimental sensitivities. Moreover, away 
from the finely tuned condition, the small masses of the light neutrinos enforce a dramatic suppression of the mixing 
(a typical seesaw behavior) and the probabilities drop by several orders of magnitude to unreachable levels. 

Model II, in contrast, requires no finely tuned condition to reach its maximal amplitudes. In this case, because 
of lepton number conservation, the heavy Majorana neutrinos come in degenerate pairs, thus acting effectively as 
Dirac neutrinos. As a result, the Majorana boxes cancel in the amplitudes and the standard Dirac boxes, which are 
generally smaller, are the only contribution. Therefore, the cancellation of the Majorana boxes occurs only in model 
II and it is due to interference between the different neutrino flavors (i.e. the flavor mixing matrix), and not because 
of Fermi statistics, as was stated in ref. 0. Despite the cancellation of the Majorana boxes, Model II gives the largest 
possible conversion probability, because the mixings are free from the seesaw suppression of model I imposed by the 
light neutrino masses. Overall, the effective coupling G mi ^/Gf in model II can reach 10~ 5 — > 10~ 4 for natural values 
of the parameters and neutrino masses around 10 TeV. 

We also included the calculation of free /i + e~ — > ~e + scattering at high energies, because this process tests the same 
amplitude as muonium-antimuonium conversion. However the cross section is so small that muonium-antimuonium 
conversion is still a more hopeful alternative from the experimental point of view. 
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APPENDIX A: THE MODELS 



In this Appendix we will briefly describe two distinct seesaw-type models, which are used in our numerical evaluation 
of the muonium-antimuonium conversion probability. Basically, these models have extended neutrino sectors as 
compared with the SM. The entire lepton sector has the same spectrum as the Standard Model (we consider Nl 
families in general) plus Nh extra right-handed neutrinos (j — 1, . . . , Nh)- The lepton sector is thus divided in 
Nl SU(2) L -doublets L (i \ N L charged singlets 1$ and N H neutral singlets 



L(i) = [% . tfP > (i = l,-..,N L ; j = l,...,N H ) . (Al) 




Concerning the interactions, only those involving charged currents are relevant to our case (lepton currents coupled 
to Ws and, in a general gauge, to charged Goldstone bosons G ): 

C lnW {x) = (~^=) [T(xh^ L (x)] W-(xr +h.c. , (A2) 

Ci nG {x) = -jJjj-^j G~ (x) l(x)m\v L (x) - T(x)m D v R (x) + h.c. . (A3) 

In Eqs. (|A2|) - (|A3|) . we used the following simplified notation: I is a column with the Nl charged leptons, i, = 1$ +1$; 

i>l is a column with the Nl standard neutrinos, v L , ^R is a column with the Njj extra neutrinos Vr\ mi is the 
standard Nl x Nl mass matrix of the charged leptons (which can be taken to be diagonal and positive, without loss 
of generality) ; nip is a new Nl X Nh mass matrix of Dirac type for the neutrinos, g is the SU(2)l coupling constant 
(g 2 = 8G F M^/V2). 

Expression (|A3I) is obtained from the Yukawa terms for the leptons and the Higgs doublet <!> (or its conjugate 
$ c = ir 2 $*): 

N L N H N l 

C Y (x) = ^LW^^WR^^^W+^^LW^^W^^W+Lc. (A4) 

i,j=l 3=1 i=l 

The first term is standard and gives the mass to the charged leptons in terms of the Yukawa couplings when 
the Higgs acquires a vev ($) = (Q,v) T /V2, namely [ m i]ij — T^ijv/^/^': it also generates the first interaction term 
in Eq. (|A3(1 . The second term is extra physics; in a similar way it generates Dirac type masses for the neutrinos, 
[mojij = Liijv/V^ and the second interaction term in Eq. (|A3() . 

In addition, in these seesaw models there are neutrino mass terms of Majorana type, oc v£vr + h.c, which can be 
gathered in a matrix uim of dimension Nh X Nh- Consequently, the totality of neutrino mass terms can be expressed 
in the form: 



-C 



\ (vz,l%) M [ < ) + h.c. = \ {vz^t) ( ° T mD )( V S ) + h.c. . (A5) 



7 1 

The superscript "c" denotes charge-conjugated fields ip c = Ctp , where C — ~ij 2r y° in the Dirac representation. The 
block in the upper left corner is zero due to the p-parameter constraints (see e.g. Ref. |l4j). The total mass matrix 
M. appearing in Eq. i|A5(l has dimension (Nl+Nh) x (Nl+Nh) and is a symmetric, in general complex, matrix with 
a texture leading to a seesaw mechanism upon diagonalization. This matrix can always be diagonalized by means of 
a congruent transformation involving a unitary matrix U (IE) 

UMU T A* = M d ■ (A6) 

Here, Md is a nonnegative diagonal matrix, and A* is a diagonal unitary matrix: (A*)y = SijX*, where A.; are complex 
phase factors (|A,*| = 1) such that A4d is made real and nonnegative. The diagonal A* matrix in Eq. (|A6|) can be 
incorporated into U by the redefinition: U ncw — (A*) 1 / 2 ?/, where U now is again unitary; A* thus reflects a freedom of 
choice for U. The Nl + Nh mass eigenstates are Majorana neutrinos, related to the interaction eigenstates v a by 
the matrix U of Eq. (|A"BT) 

/ \ Nl+Nh 

[l L cj = E u ia n ^ => f 1 i = >: r - k »«>■ • (AT1 
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The first Nl mass eigenstates are the light standard partners of the charged leptons, while the other Nh eigenstates 
are heavy (if the seesaw mechanism takes place). The factor (A*) a — A* is now recognized as the creation phase factor 
0,01 of the Majorana neutrino n, (= niL+nw), in the sense that n\ = \*n i: and more specifically {nu,)° — \*nm- 
The NlX (Nl+Nh) -dimensional mixing matrix B for charged current interactions is defined as 

B M = £/* = ([/% . (A8) 

When CP is conserved, the coefficients Bu can be chosen real and all the phase factors A^ real as well 0|: A^ — 
fjcp(jii)/i = ±1, where fjcp(jii) — ±i is the intrinsic CP parity of the Majorana neutrino m. If the symmetric mass 
matrix M. is real, then there exists a unitary real (i.e., orthogonal) matrix U — O such that 

OMO T = V , (A9) 

where T> is a diagonal and real matrix; in general, some of its diagonal elements are negative, but the absolute values 
are the same as those of Ma = T>A*: {Md)u = T>a\* — \T>u\. We see that A$ is here ±1 (= sgnpD^]). Therefore, when 
the symmetric mass matrix M is real, we are in a situation with no CP violation: Bu = On are real, and Ai = ±1. 
Even more, in that case Oij for j > Nl are also real. Thus, the reality of M. is sufficient for having CP conservation. 1 
We will see later that in this case the values of A^ either as +1 or —1, separately for each i (i = 1, . . . , Nl + Nh), 
influence our results. 

Inserting the transformations (|A7|) in the interaction terms <|A2fl and (|A3|) , taking also into account the mass matrix 
form l|A5(l and the transformation l)A6|l into Eq. l|A3jl , we obtain the explicit interactions of the charged bosons with 
the leptons in their mass basis: 

/ v iVt N L +N H 

C lnW (x) = ( — %)Y1 [B ij l i (x) 7li P L nj(x)]W-(x)^ + h. c . , (A10) 

V v2y i=1 ]=1 

/ \ Nl Nl+Nh 

C lnG (x) = G-(x) ]T J2 l B iM*) i^ p L - MjPr) nj (x)] + h.c. , (All) 

\ V W / 

where rrii > is the mass of the charged lepton Ij. and Mj > the mass of the Majorana neutrino nj. Expressions 
(|AT0|) and JATT|| were obtained already in Ref. and used in Refs. 0,^2 an d others. 
In this work, we have used two specific models: 

Model I: This is the usual seesaw model described above, with Nl left-handed neutrinos i/{l and an equal number 
of right-handed neutrinos I'm: Nh = Nl- The neutrino mass terms are those given in Eq. 1|A5(1 . The square (Nl~xNl) 
matrices rriM and tod are supposed to be invertible. 

The heavy-to-light neutrino mixings (s L 1 ) 2 = ^2 h \Um\ 2 are of the order of the squared ratio between the Dirac 
mass (mo) and the Majorana mass (ttim) scales (~ |m£)| 2 /|mj\/| 2 ). The eigenmasses of the light neutrinos are 
m vii ght ~ nT' 2 j/ rn M- Severe constraints on the |to.d| <C \itim\ hierarchy in Model I are imposed by the very low 

experimental bounds m Vu ht ~ 1 eV. Another set of constraints on the model is imposed by the present experimental 

bounds on the heavy-to-light mixing parameters (s L 1 ) 2 ~ |m._D| 2 /|mAf| 2 (~ 10 -2 , see below). 

Model II: This model is formally again of the seesaw type, because the mass matrix M. has a form similar to \hh\ . 
However, now Nh = 2Nl- This model contains an equal number Nl of left-handed (S^l) and right-handed (vm) 
neutral singlets 0,0], and thus the right-handed column Vr with Nh components is now (vr,S l c ) t . The form 
of the mass matrix M. in this model ensures conservation of total lepton number, but lepton flavor mixing is still 
possible. The neutrino mass terms, after electroweak symmetry breaking, have the form 



1 




^ ss = ^[^r,Sl)M y R |+h.c. , M=\ m T D m T M \ , (A12) 



m D 

T n _T 



m M 



which is formally of the form (|A5(I . but with a specific Nl xNl texture: the Dirac mass submatrix now has a specific 
Nl*2Nl form (m£>,0), where mo is a square (NlxNl) matrix, and is the NlxNl zero matrix; and the Majorana 



1 It is not clear to us what are the necessary and sufficient conditions on the elements of Ai for having CP conservation. If only the upper 
left Nl x Nl block of M is real, it appears that we can have CP violation, i.e., the elements Bu = U* t (I = 1, . . . , Nl) apparently 
cannot always be chosen real in such a case. 
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mass submatrix mjvf of (|A5|) is now replaced by a 2Nl x 2Nl matrix in the lower right part of M. QA12JI ■ As given 
in i|A12[l . the matrix M. gives for each of the Nl generations a massless Weyl neutrino and two degenerate Majorana 

neutrinos [2£j, Consequently, the seesaw-type restriction m Vlight ^m D /m,M ~ 1 eV of Model I is absent in Model 
II in its unperturbed form l|A12|) . Nevertheless, the experimental bounds on the heavy-to- light mixing parameters 

( s i) 2 ~ |wd| 2 /|tom| 2 (~ 10~ 2 ) do impose a hierarchy \mn\ < |wim| between the Dirac and Majorana mass sectors 
and, in this sense, the model remains of a seesaw type, but the hierarchy \mo\ < \tum\ here is in general much 
weaker than in Model I. Interestingly, nonzero masses for the Nl light neutrinos can easily be included in Model II by 
introducing small Majorana mass terms for the neutral singlets Sil 1 i.e., small nonzero elements in the NlxNl lower 
right block of M. The mixings of heavy-to-light neutrinos are not significantly affected by these small corrections. 

The above two models, for specific choices of the Nl x Nl matrices mo and vtim , give us specific values of the 
mixing-matrix elements (|A8() and eigenmasses Mj, i.e., the parameters which then determine the charge current 
interactions of IjAlOp and (|A11(1 . 



APPENDIX B: BOX DIAGRAM AMPLITUDES WITH MAJORANA NEUTRINOS 

Here we outline the derivation of the expression for G M ^ in Eq. (|10[l at one loop order (box diagrams) , defined by 
the effective lagrangian of Eq. 

There are basically four types of box diagrams: two containing generic neutrinos (Fig. 1) and two exclusive for 
Majorana neutrinos only (Fig. 2). According to the contractions of the fields with the external states, the diagrams 
of Fig. l.b, 2. a and 2.b have the same relative sign, while the diagram of Fig. l.a has the opposite sign. 

There are also basically two momentum integrals encountered in the box diagrams, in the limit where external 
masses and momenta are neglected compared to the loop momenta. The first one is: 

Iij = / (2^ 4 (q*-M?)(q*-Mj)(q*-M*,)* = (^M 4 /^'^' (B1) 
where we denote Xj = Mj/M^- and where Jixi^Xj) is a dimensionless expression: 

1 { ( 1 a;,- In Xi 



J ^ = -(^H(i-^) + (i- ~ >• (B2) 



The second integral is: 



d 4 q q 2 

(27T 

where the dimensionless expression is 



K " - I (2„)« (^-M?)(,'-M?)( 9 =-M5,)= {infMl^" 1 '^ (B3) 



" " { (ih; + - - "<>} ■ < B4 > 

The transition amplitudes T^fy and T^^, which correspond to the diagrams of Figs, l.a and l.b when both 
internal bosons are W, are: 

' T wl \ = (-4*} B* s Bt 4 B.,iB 



x f _A 1 f ( + 1) [Tiik'h^uPLvip')} [v{k)Yi^PLu{p)] \ , , 

J (2tt)4 (q^-M^-Mf^-M^y \ (-1) [u(k') 7 J lv P L u(p)} [v(k)l l 'ft* S ^W)] j K ' 

Overall we label p and p' as the momenta of the initial and final electron (positron) , and k and k' as the initial and 
final momenta of the (anti-) muons. To simplify the notation, we have omitted the spin and flavor labels in the spinors, 
since the momenta specify them unambiguously, (c.f. Fig. 1). We notice a factor (—1) in the second amplitude; its 
origin is the anticommutation of the fermion fields involved in the Wick contractions with the external states. Now, 
using the identity 

7p la Ju = it\^avl X lb + ~i^g av + Jug^a ~ laQ^w (B6) 
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the spinor structure can be simplified to the following form (we have omitted the spinors in the brackets): 

h^P L ] [i u i^p L ] = 4<?V [ la p L ] b P P L ] ■ (B7) 

Finally, after replacing q a q@ — ► g a @q 2 /A due to the Lorentz invariance of the integrals and doing the following Fierz 
transformation in the second amplitude: 

[ui la P L u 2 ] [v 3 <y a P L v 4 ] = - [u lla P L Vi] [v 3l a P L u 2 ] (B8) 

both amplitudes in Eq. I|B5[1 become identical: 

T ww = r ww = i^Ki^B^MxuXj) [u{k') la P L v{p')] [v(k)l a PLu(p)} . (B9) 

(la) 
WG ' 



In the gauge we are using, one must also include the diagrams with one and two Goldstone bosons, i.e. % 



Tqw ^ Tqq\ etc. The calculation goes through similar steps as those already stated for the box with two W's, but 



contains different propagators and therefore different combinations of the integrals IjBlfl and l|B3(l 
The sum of all these diagrams gives the total amplitude of the "Dirac" box diagrams of Fig. 1: 



T(i) - \ " IVO) , , n-{ 

1 — I WW GG G 



-(Ik) 

L GW+WG 

k— a, 6 

,2 



= - l ^k 2B* ei B* ej B^B M F Box ( Xi , Xj ) [u{k') la P L v{ P ')] [v(k)rPLu(p)} , (BIO) 

where FBox(xi, Xj) is given in Eq. and is identical to the expression Fft ox (xi, Xj) of Ref. |12| . 

The diagrams of Fig. 2 are, on the other hand, different in the sense that they exist only if the internal neutrinos are 
of Majorana type. The Wick contractions of the Majorana fields can be performed in the usual way if we previously 
transform two of the four vertex operators as follows: (we chose the vertices with the initial anti-muon and with the 
final anti-electron): 

^nixhuPLrijix) = -Xjnj(x)j u P R ip^(x), (Blla) 
ni{x)^^PL^e{x) = -Xlip^(x)j,j,P R rii(x). (Bllb) 

Here, Xj is the creation phase factor of the Majorana neutrino rij, i.e. n'j = X*rij, and the superscript "c" denotes 
the conjugated field ip c — , where C — —ij 2 ~f° in the Dirac representation. After performing the transformations 



(|Bllf) . we can contract the n(xk)'s and n(x n )'s into normal neutrino propagators (see Ref. |14|). Let us first consider 
the two diagrams of Fig. 2 with two internal W's. These diagrams differ from each other only in the way the Ws are 
contracted. The amplitudes are: 



T (2a) 

1 ww 



3 I = ("^) (^(^W^-MiC-l) 

d A g 1 [ [u{k') llllv P R u c {k)] & c (p')^YPlu{ P )] 

(2tt)4 ( g 2_ M 2 )(g2 _ M 2 )(g2 _ M 2 /) 2 | [ ¥ (fc') 7M> P flUc (fc)] [u c (p') 7 ^P LU {p)] 



(B12) 



We denoted as u c the u-spinors associated to the positively charged fermions (usually considered "antiparticles"), but 
operationally these are normal u-spinors. The factor (—1) before the integral is important and comes from the Wick 
contractions of the fermion fields with the external states, given the sign convention of the previous diagrams. Using 
the identities 

- <T - , ^75 = -e^cfiv"* , (B13) 
the spinor structure in both amplitudes of Eq. (|B12|) can be simplified to the same form: 

[WvPr] WYPl] = [lulM [Yi^Pl] = 4 [Pr] [Pl] ■ (B14) 
This is the reason why the diagrams in Fig. 2 do not cancel but, on the contrary, they add to twice their value. 
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Using the following Fierz identity: 

2 \uiPru 2 ] [usPlua] = [ui^oPlua] [m37 q -PrM2] (B15) 
we can transform the spinor product into a product of currents. The amplitudes of Eq. l|B12fl then result in: 



?ww = ?ww = -<2§r ( S -) 2 {B M f ^xjjixuxj) [u(k') la P L u(p)\ [u^p'^Pru^)] . (B16) 

Now, the diagrams with Goldstone bosons instead of W's contain scalar currents which, in analogy to Eqs. HBllfl . 
should be transformed before doing the Wick contractions according to: 

^„(x)P l ,r nj(x) = X j n j (x)P L , R yb^x), (B17a) 
n l {x)P L ^ R ip e (x) = \*ip%(x)P Lt R rii(x). (B17b) 

Besides this detail, the calculation of the diagrams of Fig. 2 with one or two Goldstone bosons replacing the W's 
follows the same steps as above. 

From Eq. (|B16|) and the analogous expressions corresponding to the diagrams with Goldstone bosons, it is clear 
that both diagrams of Fig. 2 give identical results. 

We should now transform the it-spinors of the positively charged fermions into v-spinors, as it is customary, using 
the following conjugated current identity: 

u c (p')7 Q ^Wfc) = v(k)j a P L v(p') (B18) 



and use the Fierz identity (|B8|) to arrive at the same spinor structure as in the amplitude of Fig. 1, namely Eq. 1|B10|I 
The total result for Fig. 2 is then: 

/t-(2) _ \ |V(kO , /t-(«) , <7-(«) 
1 — I WW + 2 GG 1 GW+WG 

K—a.b 
2 

' (B* ei ) 2 (B M ) 2 X 3 X* G Bo ^x llX] )[u(k'h a P L u(p)][mrPLv(p')] , (B19) 



where the function Gbox is given in Eq. (j!2(l . which is identical to that of Ref. (l2^ . 

The sum of all diagrams is r(!) + T@>[ Eqs. l|B10)l and from this amplitude one deduces the efective 

coupling G M jtj of Eq. (fT*0|l . 
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